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Abstract 

The effects of supersymmetry breaking are usually parameterized by soft couplings 
of positive mass dimensions. However, realistic models also predict the existence of 
suppressed, but non- vanishing, dimensionless supersymmetry-breaking couplings. These 
couplings are technically hard, but do not lead to disastrous quadratic divergences in 
scalar masses, and may be crucial for understanding low-energy physics. In particular, 
analytic scalar quartic couplings that break supersymmetry can lead to intermediate scale 
vacuum expectation values along nearly-flat directions. I study the one-loop effective 
potential for flat directions in the presence of dimensionless supersymmetry-breaking 
terms, and discuss the corresponding renormalization group equations. I discuss two 
applications: a minimal model of automatic .R-parity conservation, and an extension of 
this minimal model which provides a solution to the /i problem and an invisible axion. 



1 Introduction 



The primary motivation for supersymmetry (for reviews, see [jl], H|) as an extension of the 
Standard Model is that it can stabilize the hierarchy associated with the electroweak scale. This 
implies that supersymmetry is softly broken [||, ||. The technical version of this requirement is 
often taken to mean that quadratic divergences in radiative corrections to scalar masses must be 
absent to all orders in perturbation theory. However, the determination of which couplings are 
soft then depends on whether the theory contains any gauge-singlet chiral superfields which can 
engender tadpole loop diagrams. This is a rather obscure criterion from the low-energy point of 
view, since the minimal supersymmetric standard model (MSSM) contains no such fields, but 
reasonable extensions of it often do. Furthermore, if supersymmetry is spontaneously broken, 
then the low-energy theory will generally include couplings that are not soft according to the 
technical definition. In realistic models, the hard supersymmetry breaking couplings are usually 
expected to be highly suppressed, and that is why they are traditionally neglected. In this paper, 
I will explore in detail some circumstances in which "hard" supersymmetry breaking couplings 
nevertheless are important, and even crucial, for understanding physics at low energies. 

Let us suppose that the spontaneous breaking of supersymmetry can be parameterized by 
the vacuum expectation value (VEV) of an auxiliary F-term for some chiral superfield X. The 
appearance of supersymmetry breaking in the low-energy theory can then be understood as 
coming from non-renormalizable Lagrangian terms which couple X to the other chiral super- 
fields <I> and gauge field-strength superfields W£ in the theory. These terms are suppressed by 
powers of some large mass scale M, which in supergravity-mediated supersymmetry breaking is 
related to the Planck mass. The relevant couplings of lowest dimension include, schematically^ 

- c = + yi [x ^ ]f + Ti [x ® 2]f ) + cx - + w [x * x ^^ ]d - (L1) 

So, with (X) = 99F, one finds the usual soft terms 

- £ = (s A ° A ° + s 03 + f* 2 ) + c ' c ' + iJ^ (L2) 

corresponding to gaugino masses, cubic scalar couplings, analytic scalar squared masses and 

non-analytic squared masses. We can therefore make the identification F/M ~ mw, indicating 

that these soft terms are very roughly of order the electroweak scale. The analytic scalar 

squared mass is parameterized here as being proportional to a corresponding superpotential 

mass parameter /x. This is not strictly necessary, but corresponds to the phenomenologically 

^The subscripts D and F in eq. ( [L.l| ) correspond to integrations J d 4 # and J d 2 respectively. Indices a,b,c . . . 
are used for gauge adjoint representation indices, a — 1, 2 is a two-component fermion index, and indices i,j,k... 
will be used for chiral supermultiplet labels. The chiral superfields $ have complex scalar and fermion components 
4> and tp, and X a denote gaugino fields. 
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viable picture of the MSSM in which \i is itself expected to be of order F/M or 1 TeV, much 
less than the high scale M. (In a more complete model, [i is likely to be replaced by the VEVs 
of some other fields so that it is related to F/M in some way.) 

However, there is no good reason why the only terms that couple X to the other fields 



should be the ones shown in eq. (L_l). All possible renormalizable supersymmetry-breaking 
terms can arise from appropriate non-renormalizable supersymmetric couplings involving one 
or two powers of X. This is shown in Table 1, in which the interactions are written schemat- 
ically along with their possible origin and softness according to the technical definition. I list 
each coupling as either "soft" (if it never leads to quadratic divergences in conjunction with 
renormalizable supersymmetric couplings), or "maybe soft" (if no quadratic divergences can 
be induced provided that there are no gauge-singlet chiral superfields), or "hard". The last 
column lists the lowest-dimension possible origin for the term. The subscript D or F indicates 
whether the original term involves / d 4 9 or J d 2 9. The penultimate column then lists the re- 
sulting naive suppression, assuming that the coefficient of the original term involving X, X* is 
of order unity. Of course, this last assumption can easily be violated by symmetries, so that 
the actual suppression may be stronger or weaker than the naive estimate shown. 

The existence and potential importance of the non-analytic cubic couplings 4> 2 4>* and chiral 
fermion mass terms ifiip have been recognized in several papers e.g. H [5|, |6|, [?], ||, [|. These 
terms are closely related to each other, since any chiral fermion mass term ipip can be absorbed 
into the superpotential, at the cost of also redefining the </> 2 c/>* couplings if there are superpo- 
tential Yukawa couplings present. In this sense, hard chiral fermion mass terms are redundant. 
The chiral fermion gaugino mass mixing terms Xip can exist if a chiral fermion is in an irre- 
ducible representation which is also found in the adjoint representation of the gauge group ||. 
This does not occur in the MSSM, but can easily happen in extended models. The one-loop 
renormalization group (RG) equations for all of these dimensionful supersymmetry breaking 
couplings have been recently worked out in ref. ||. 

The dimensionless hard supersymmetry breaking couplings include three distinct types of 

quartic scalar couplings (j) 4 , (fi 3 4>*, 4> 2 4>* 2 , and various types of scalar-fermion- fermion couplings 

as allowed by gauge symmetries. For example, there are ^>AA couplings if a chiral supermultiplet 

transforms as a representation found in the symmetric product of the adjoint with itself. There 

are non- analytic Yukawa couplings (jftptjj, as well as non-symmetric analytic Yukawa couplings 

<t>iil)jipk- Unlike Yukawa couplings following from the superpotential, the latter need not be 

symmetric under interchanges of i,j or of i, k. Note that all of the couplings listed in Table 1 

that have a D-term origin are suppressed byQ \F\ 2 /M 4 , while those that have an F-term origin 

'If D-term VEV(s) plays an important role in supersymmetry breaking, then \F\ 2 /M 2 can be replaced (or 
added to) by D everywhere in Table 1. 
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Table 1: Classification of all renormalizable super symmetry breaking interactions. Chiral sca- 
lars and fermions are represented by and ip, and gauginos by A. The last column indicates the 
lowest dimension operator which can give rise to the term through spontaneous supersymmetry 
breaking with (X) = OOF and a high scale of suppression M. The resulting naive suppression 
is shown in the third column. 



Type 


Term 


Naive Suppression 


Origin 




#* 


\F\ 2 2 




soft 




^ ~ iJ,m w 






<P 3 


jj ~ m w 


^[*<& 3 ]f 




XX 


jt ~ m w 








IF M 




maybe soft 




IF M 








JFP _ rn^ 
IF M 


-± s [XX*D a $W a ] D 




4 


IF M 






3 0* 


|F| 2 m 2 , 
IF IF 


^[XX*^*] D 




2 0* 2 


|F| 2 m 2 , 

IF IF 


-^[XX*$ 2 <S>* 2 ] D 




0-0-0 


IF IF 


^[XX* $D a $D a $] D 


hard 


0*-0-0 


IF IF 


-±z[XX* $*D a <f>D a <f>] D 




0V>A 


IF IF 


j^[XX* *D"*W a \ D 




0*^A 


IF IF 


j^[XX* $*D a <f>W a ] D 




0AA 


IF M 


jfa[X$W a W a ] F 




0*AA 


IF IF 


^[XX* <Z>*W a W a ] D 
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are only suppressed by F/M 2 . 

There are at least two reasons why the suppressions of the "maybe soft" and "hard" cou- 
plings listed in Table 1 need not render them irrelevant. First, it may be that our simplest 
notions of spontaneous supersymmetry breaking are incorrect, so that the naive suppressions 
that are usually implicitly assumed are not present. For example, it may be that the high scale 
M which governs the suppressions is actually far below the Planck scale. Or, several different 
high scales, some much smaller than others, may govern the suppression of these terms. Second, 
the MSSM and other realistic supersymmetric models generically have many "flat directions" 
[10, [TlJ 12, [HJ along which the renormalizable supersymmetric part of the scalar potential 



vanishes identically. Along the fiat directions, the effects of the usual dimensionful soft terms 
become less relevant at intermediate and high scales, so that the dimensionless supersymmetry- 
breaking couplings become significant despite their suppression. 

In this paper, I will concentrate on the dimensionless supersymmetry breaking couplings of 
the type </> 4 . While these couplings are technically hard,| they are soft in the practical sense 
that they do arise generically in models with spontaneous supersymmetry breaking. They 
are typically suppressed only by one power of F/M, unlike the cj^cff and cj) 2 (f)* 2 dimensionless 
terms and the 4> 2 4>* scalar cubic term. Section || contains the one-loop RG running for these 
couplings, and gives analytic solutions for the case that gauge couplings dominate, or more 
generally when the </> 4 term parameterizes a flat direction of the renormalizable supersymmetric 
lagrangian. These results are useful for relating high-scale boundary conditions on the soft terms 
to the possibility of symmetry breaking at much lower scales. In section |3L I will describe how 
the supersymmetry-breaking 4 couplings are relevant for symmetry breaking at intermediate 
scales, including a discussion of the one-loop effective potential. In sections ||] and |j I will apply 
these considerations to two examples: a minimal model of B — L breaking to i?-parity, and an 
extension of this minimal model which also incorporates an invisible axion and a solution to 
the n problem. Section ^ contains some concluding remarks. 

2 Running of supersymmetry-breaking 4 couplings 

In this section I will discuss the RG equations relevant for running an /i^ 4 -type coupling from 
a high scale (where a boundary condition on it is to be provided) to an intermediate scale or 
other scale of interest. Consider a superpotential of the form 

w = -Lx^i^SfcS, + iy^i^* , (2.i) 



§ Actually, <jf> 4 couplings are always technically soft at one- loop order; they only give rise to quadratic diver- 
gences in scalar squared masses at two-loop order or higher. 



5 



with a gauge group with coupling(s) g a . There are also super symmetry breaking couplings 

" ^breaking = {^^MjMl + ^'Vi^fc + ^a^X^ + C.C. + (m 2 )&*% (2.2) 

where M a is the gaugino mass, (m 2 )l is the scalar squared mass, a l ^ k is the holomorphic soft 
(scalar) 3 coupling, and h l ^ kl is the super symmetry-breaking (scalar) 4 coupling. Note that both 
the supersymmetric coupling x and the corresponding supersymmetry-breaking coupling h are 
defined with a 1/M factor where M is the high mass scale. Therefore x is dimensionless, and 
h has dimensions of mass and should be roughly of order mw- The RG equations can be 
computed by requiring that large logarithms involving the cutoff can be absorbed into coupling 
constant and field redefinitions, introducing a renormalization scale Q. The /3-function for any 
running parameter is equal to its derivative with respect to t =ln(Q/Qo) where Qq represents 
some fixed energy scale. 

As a point of reference, consider first the well-known one-loop beta functions of the renormal- 
izable supersymmetric and soft couplings. For the gauge coupling and gaugino mass parameter, 
one has 

16tt 2 /3( 5 2 ) = ba gl- 16vr 2 /?(M a ) = 2b a g 2 a M a , (2.3) 

with 

b a = S a (R) - 3C(G a ), (2.4) 

where S a (R) is the Dynkin index summed over all chiral supermultiplets and C(G a ) is the 
Casimir invariant of the adjoint representation. The normalization is such that for SU(N), 
C(G a ) = N and a fundamental representation contributes 1/2 to S a (R). The superpotential 
Yukawa interactions obey 

1 



16tt 2 /3(Y^ 



Y^Y pqn Y^ k - 2g 2 a C a (i)Y^ 



+ + (i<-fc), (2.5) 



and the scalar squared masses have beta functions 



Wn 2 f3({m 2 )l) = ^Y tpq Y^ n (m 2 y n + ^Y^Y pqn (m 2 )2 + 2Y ipq Y^{m 2 f r + a ipq a^ 

-85{g 2 a C a (i)\Ma\ 2 + 2 ! fa?Tr[t a m 2 \ . (2.6) 

where the last term explicitly involving the generators of the gauge group t a vanishes for non- 
abelian groups. Finally, the beta functions for a l ^ k terms are given by 



167r 2 /3(a 



ijk\ _ 



Y l3n Y nvq aP qk + \Y^Y pqn a^ k - 2g 2 a C a (i)a^ k + ±g 2 a C a (i)M a Y^ k 
+ + (i++k). (2.7) 
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Figure 1: The terms in the RG equation fl2.8| ) proportional to gaugino mass come from this 
logarithmically divergent Feynman diagram. The dashed lines correspond to scalar fields, the 
solid lines to chiral fermions, and the solid line with a wavy line superimposed corresponds to 
a gaugino line with a mass insertion. 

These results employ the standard convention that Yijk = (Y* 3 )*. If there are several distinct 
gauge couplings g a , then a sum over the index a is implicit where appropriate. 

Now with the conventions outlined and illustrated above, consider the RG running of the h 
and x couplings. The coupling h i: ' ' runs with renormalization scale according to 



lQit 2 (3{h ijkl ) 



1 1 1 

ZhVPIV ~v/~Tikl I Uikpqyr ynjl j UilpiV ynjk 

2 " 1 pqn 1 2 ^ n 2 



+-Y^Y pqn h^ kl - 2g 2 a C a (i)h^ kl + 4g 2 a C a (i)M a x^ kl 



2 

+{i j) + (i <h> k) + (i <->■ I). (2i 



The first three terms in eq. ( |2.q ) are due to vertex renormalization, indicative of the fact 
that supersymmetry has been broken. The next two terms have the standard superfield wave- 
function renormalization form. However, the last term is somewhat unusual in that it involves 
a non-renormalizable coupling (x l ^ kl ) contributing to the RG running of a renormalizable cou- 
pling. This term comes from the logarithmically divergent Feynman diagram shown in Figure 
[l|. Usually, such contributions are simply ignored, because the non-renormalizable coupling is 
understood to be suppressed by the high scale M. However, in the present case it cannot be ne- 
glected because h is also expected to be suppressed. It is easy to see that it is of the same order 
as the other terms, given the estimates h/x ~ F/M ~ M a ~ mw- We will see further evidence 
for the necessity of such terms, and the consistency of neglecting other non-renormalizable 
effects, in section p| 

The RG running of the corresponding non-renormalizable superpotential parameter x l ^ kl is 
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just given by superfield wavefunction renormalization: 

+ (i<->j) + (i*->k) + (i<->l). (2.1 



16ir 2 p(x ijkl ) = -Y ipq Y pqn x njkl - 2g 2 a C a (i)x ijkl 



Now, typically the he/) 4 terms are important when investigating directions in field space 
in which the supersymmetric part of the scalar potential is both D-flat and F-flat at the 
renormalizable level. For this purpose, we can restrict our attention to cases in which the 
Yukawa couplings Y l i k do not connect any two fields involved in the flat direction. In general, 
D-flat directions correspond to analytic polynomials in the (pi, which in this case correspond to 
non-vanishing terms h l i kl ^(pjCpk&l i n the scalar potential. This means that when investigating 
such flat directions, the first three terms of eq. ( |2.8| ) are absent. In many cases, the wavefunction 
renormalization factors Y ipq Y pqn will also be proportional to 5 l n , so that each component of h^ kl 
as well as x l ^ kl will evolve separately under RG running. (An obvious special case of this occurs 
if the gauge interactions dominate over the pertinent Yukawa couplings.) If so, then the ratio 
of each component of an coupling to the corresponding x^ coupling will obey a simple 
RG equation: 

J t Wx) = ^Y,9lC t (a)M a + (i - j) + (i -> k) + (i - /), (2.10) 

a 

where we omit the indices i, j, k, I in the ratio h/x for simplicity. The point is that the superfield 
wavefunction renormalization factors cancel out of the RG running for the ratios h/x, leaving 
only the terms proportional to gaugino mass. This is particularly useful since, as we shall 
discuss in the next session, the ratio h/x is most important in deciding whether a non-trivial 
minimum can occur along a flat direction. Equation (2.10) is easily solved in conjunction with 
eq. fl2.3j ), with the result 

h/x = ho/x - Yl <9l fi M afi C a (i) + (i - j) + (j _ fc) + (i _ J), (2.11) 



where 



-1/2 

K a (t) 



l + (t -t)^b a 



(2.12) 



and g a fi, M a ^, and ho/xo are the values of g a , M a , and h/x at the reference scale to- 

For example, one might choose to specialize further and assume that the gauge couplings 
and gaugino masses unify to values gu and m\/2 at a unification scale to. Then at lower scales 
one has the results: 

9a = K a gu ; M a = K 2 a m 1/2 , (2.13) 
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If the soft scalar masses have an initial value m\ at the scale %, then at a lower scale t they 
are given by the well-known result 



m 2 



ml + 2m\ /2 £(1 - 4)C a (i)/b a + ..., (2.14) 



where the ellipses denote negative-definite contributions from Yukawa couplings. The ratio h/x 
at the scale t is just 

h/x = h /x - (t - i)0m 1/2 + (* ~> ■?) + (* *0 + (* ~> 0- ( 2 - 15 ) 



The result eq. ( 2.1 5| ) implies that the running of h/x can have either constructive interference, 



if ho/xo and m^ii have opposite signs in our convention, or destructive interference if ho/xo and 
m\j2 have the same sign. (In general, h/x and mi/2 can be related by a non-trivial complex 
phase, of course.) In the case of constructive interference, the dimensionless ratio h/xm\n can 
become quite large in the infrared even in the case of a flat direction made up of fields charged 
under an abelian symmetry,^ as we will see in Section 4. This is the situation that most favors 
a non-trivial VEV at an intermediate scale. 

3 Symmetry breaking along flat directions 

A toy model which illustrates the essential features we want is as follows. Consider a pair of 
chiral superfields which have opposite charges ±1 under a gauged U(l) symmetry with coupling 
g. Assuming that there is a flat direction in the renormalizable, super symmetric part of the 
scalar potential, the leading contribution to the superpotential can be written as: 

W = -^jS 2 S 2 . (3.1) 

The coupling x is taken to be real and positive without loss of generality. There is also a U(l) 
D-term contribution, so that the supersymmetric part of the scalar potential is 

Vsusy = ^\SS\\\S\ 2 + |5| 2 ) + 9 ^(\S\ 2 - |5| 2 ) 2 . (3.2) 

The absence of a superpotential mass term SS can be ascribed to a discrete symmetry, or 
perhaps viewed as a general result of superstring models which generically do not have tree- 
level masses. After spontaneous supersymmetry breaking, there are also contributions to the 
scalar potential: 



VI 



breaking 



m 2 |5| 2 + m 2 |5| 2 +(A^5 2 + c.c). (3.3) 



^The growth of h/xm^/2 can be even more dramatic if the fields participating in the flat direction are charged 
under an asymptotically free gauge interaction; see Appendix A for an MSSM example. 
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For small field strengths s, the soft scalar masses dominate. If m 2 + m 2 > 0, then s = is a 
local minimum of the potential. However, the /i-term contribution to the potential can always 
be made negative by a suitable rephasing of the fields, without affecting the phases of any of 
the other terms. Therefore, the presence of the h couplings always favors the existence of a 
non-trivial minimum for S and S, with spontaneous breakdown of the U(l) symmetry at an 



intermediate scale of order yJmy/M along the D-flat direction. For the largest field strengths, 
the superpotential terms dominate and stabilize the minimum, ensuring that the potential is 
bounded from below. 

The .D-term contribution in eq. ( |3.2[) just forces S and S to have VEVs of very nearly equal 
magnitude. Therefore, we can first evaluate the scalar potential in the approximation that the 
minimum occurs exactly along the flat direction \{S)\ = \{S)\. (Later in this section we will 
discuss the important issue of deviations from D-flatness.) Parameterizing the flat direction by 
s 2 = \S\ 2 = | S | 2 and S 2 S 2 = s 4 e ie , the tree- level potential takes the form 

1 2r 2 
V (s,9) = (m 2 + m 2 )s 2 + —(he w + /^V + j^s e . (3.4) 

A non-trivial local minimum, if one exists, occurs for = 9 m [ n with 

he* 9 ™ 3 * = h*e~ i0 * i * = -\h\. (3.5) 

Using the estimates h ~ mw and m 2 , m 2 ~ m 2 ^ , and x of order unity, one can check that all of 



the terms in eq. ( |3.4D are comparable, of order m^M, when s ~ ^my/M. The condition for a 
local minimum of the tree-level potential is then easily seen to be [14, ^, 15]: 

\h\ 2 - 6x 2 (m 2 + m 2 ) > 0, (3.6) 

with 

.2_„2_ M U , TJHZi ; my/M. (3.7) 



s = s 



°~ 6x 2 



H + V l^l 2 - 6x 2 (m 2 + m 2 



Note that for smaller values of x the location of the minimum is pushed to higher scales, but 
only the ratio \h/x\ is important in deciding whether a non-trivial local minimum exists. 

In sections |] and [5], we will discuss applications in which non-trivial minima can arise at 
intermediate scales in the way just described, modulo certain minor complications. Now, a 
complete theory of supersymmetry breaking should predict the values of the supersymmetry 
breaking parameters including m 2 , m 2 , and h/x, and therefore in principle should be able to 
tell us whether or not a non-trivial minimum actually exists. If one has in mind a supergravity- 
mediated type of model for supersymmetry breaking, these predictions will take their simplest 
form for running parameters near the Planck scale, and will have to be run down to the in- 
termediate scale. One way to understand this explicitly is to construct the one-loop effective 
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potential. The result of doing this is that h, x, m 2 , and fn 2 in the above discussion should 
be replaced by running parameters evaluated near the scale of the possible VEV, plus small 
calculable one- loop corrections. 

It seems worthwhile to see how this goes explicitly in the present example by constructing 
the effective potential along the flat direction. First, it is important to note that despite the 
presence of the hard supersymmetry-breaking coupling h, at one loop-order there is no field- 
dependent quadratically-divergent contribution to the scalar potential, since 

STr[M 2 } = 2{m 2 + m 2 ) -2\M X \ 2 , (3.8) 

even when deviations from D-flatness are allowed. Here Ai 2 are the eigenvalues of the scalar- 
field-dependent (mass) 2 matrix for the real scalars, two-component fermions, and vector bosons 
in the theory. M\ is the U(l) gaugino mass. The supertrace STY denotes a sum over these modes 
weighted by (— l) s (2s + 1) where s is the spin of the particle. Equation (|3.&j ) reflects the fact 
that analytic h(f> 4 couplings are technically soft at one-loop order. Field-dependent quadratic 
divergences do arise at two-loops, but they are of order AV ~ A 2 \h\ 2 c/)*(j)/M 2 multiplied by a 
two-loop phase-space factor and are therefore safely negligible in the following one-loop order 
calculation, as long as the ultraviolet loop momentum cutoff A is not much larger than the 
high-scale suppression mass M. 

There are subtleties involved in evaluating the one-loop effective potential for arbitrary S 
and S, since the lagrangian contains non-renormalizable couplings which lead to divergences 
in the non-quadratic part of the Kahler potential. We will return to those issues briefly at the 
end of this section. However, to a good approximation, we can reliably evaluate the effective 
potential near the flat direction including only one-loop contributions proportional to g 2 . Then 



the one-loop effective potential is given in the DR scheme [16| and in Landau gauge by the 
usual expression V = Vq + V\ , where 



Vt = gSTr {M 4 [ln(M 2 /Q 2 ) - 3/2] } , (3.£ 



64tt 2 

with Q the renormalization group scale. 

The important contributions of order g 2 to the one-loop effective potential come from loops 
involving the massive vector supermultiplet which arises when S and S are given VEVs. The 
massive vector supermultiplet fields consist of one real scalar, a pair of two-component fermions, 
and a gauge boson. They have squared masses, respectively, 

Ml = 4g 2 s 2 + l( m 2 + m 2 )-^(he^ + h*e-^s 2 -^- 2 s 4 -, (3.10) 

M\ /2 = Ag 2 s 2 ±2gs\M^-^s 2 \ + l -\M x \ 2 + ^ I - 2 s\ (3.11) 
M\ = 4g 2 s 2 , (3.12) 
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up to terms of order m^/M. Putting these into eq. ( |3.9| ) and expanding while keeping only 
terms which can contribute proportional to g 2 in the scalar potential, one finds: 



Vi(s,9) 



16vr 2 



(4|M A 



+ln 



'Ag 2 s r 

. Q 2 



(m 2 + m 2 -8\M x \ 2 )s 2 



10x 2 



+—(4xM x e ie + 4xM{e 



he 



ie 



M 2 
/i* e - i9 )s 4 



As a check, we can require that the potential Vq + V\ is RG-invariant: 



/?(*)!- + (3{h)-0r + PW-JL + /3(m 2 + m 2 ) 



a 



dx dh 
This equation is solved by 



dh* 

0(x) 
(3(h) 
I3(m 2 + m 2 ) 

Is 



d(m 2 + m 2 



d 

sls dS 



v- dVl 



2vr 2X ' 
JLj(2xM x -h); 

-4l^l 2 ; 

9 2 



16vr 2 ' 



(3.13) 

(3.14) 

(3.15) 
(3.16) 
(3.17) 
(3.18) 



In particular, eq. ( [3.16 ) confirms the RG running for the h(p^ couplings found in Section 2. 

Now one can minimize the one- loop effective potential Vq + V\ with respect to s and 9. The 
result is that a non-trivial local minimum exists provided that 



\h\ 2 -6x 2 (m 2 + m 2 



1 + 



12^ 2 



+ 



24vr 2 



\h + 6M x x\ 2 > 0, 



(3.19) 



with all running parameters evaluated self-consistently at Q = 2gs m { n , where s = s m i n is the 
resulting minimum. [Here we drop terms of order (g 2 /IQtt 2 ) 2 .] So, the region of parameter 
space in which an intermediate scale VEV is stable is actually increased (and typically only 
slightly) by the one-loop radiative corrections compared to the tree-level constraints, since the 
last correction term is positive-definite. The minimum is achieved for 9 = 9 m { n satisfying the 
same condition he 1 min = h*e min = —\h\ as for the tree-level potential, with h now equal to 
its running value at Q = 2gs m \ n , and with s m j n explicitly given by 

9 M 



6x' 2 



h'\ + \J\h'\ 2 - 6x' 2 (m 2 + m 2 )' 



(3.20) 
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where 

a 2 

(m 2 +m 2 )' = m 2 +m 2 - f-^|M A | 2 ; (3.21) 
h> ee (l + -£^)h + £^M x x; (3.22) 

x> ee + (3.23) 

with all quantities on the right side evaluated self-consistently at Q = 2gs m m. 

The one-loop corrections indicated in eq. fl3.19|) -( j3.23 ) are typically not overwhelming. Fur- 
thermore, the exact location of s m ; n is not known unless the overall magnitude of x is known, so 
without quite detailed model input it is not clear at exactly what scale to impose the condition 
for a non-trivial minimum. (Specific models of supersymmetry breaking can predict the ratio 
h/x with good accuracy more plausibly than they can predict h or x separately.) Therefore, 
we are typically justified in simply imposing the tree-level conditions for a non-trivial minimum 
of the scalar potential, using running parameters evaluated at some reasonable guess for an 
intermediate scale. That is the procedure we will follow in the following sections. 

In the preceding discussion of the one-loop effective potential, we included only g 2 con- 
tributions and only field configurations on the flat direction parameterized by s and 6. It is 
interesting to consider the more general case in which these restrictions are not made. To that 
end, consider the STr[7W 4 ] which makes a contribution to the scalar potential proportional to 
InQ 2 . It is: 

STi[M 4 } = 2(m 4 + m 4 )-2|M A | 4 + 6 5 2 (m 2 -m 2 )(|S| 2 -|S| 2 ) 

+2g 2 {m 2 + m 2 - 8|M A | 2 )(|S| 2 + |5| 2 ) + 8 5 4 (|5| 2 - \S\ 2 ) 2 

+f^(16*M A - Ah)S 2 S 2 + c.c. - ^-(|S| W + \S\ 2 \Sf ) 

+j^(\h\ 2 + 2x 2 M 2 )|5| 4 + j^(\h\ 2 + 2x 2 m 2 )|S| 4 

+^[\h\ 2 + x 2 (m 2 + m 2 )]|S| 2 |<?| 2 + ^r(\S\ 2 + \S\ 2 )S 2 S 2 + c.c. 

-^(I^l 2 " \S\ 2 )(\S\' ~ \S\ 4 ) + ^(3|S| 4 |5| 4 + \Sf\S\ 2 + |S| 2 |5| 6 ). (3.24) 

The first few terms, up to (ISI 4 ^! 2 + I^PlS"! 4 ), have the same form as terms already present 
in the tree-level scalar potential. However, the remaining terms do not. In particular, the last 
terms involving (|5| 2 - |5| 2 )(|5| 4 - |5| 4 ) and (3|5| 4 |5| 4 + |S| 6 |5| 2 + |S| 2 |S| 6 ) do not involve 
any supersymmetry breaking couplings and yet appear to have no counterpart in the tree-level 
supersymmetric scalar potential. 

The superpotential is not renormalized. Therefore, to understand the appearance of these 
terms with a InQ 2 coefficient, one must include non-quadratic Kahler potential terms which 
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are renormalized and whose couplings run logarithmically with RG scale. The Kahler potential 
can be expanded in powers of 1/M according to: 

K = S*e 2 ^ x S + S*e- 2 ° X S + -h- S * 2 e^ x S 2 + ^S*S*SS + A-S* V 4 ^S 2 + . . . (3.25) 

AM 2 M 2 AM 2 

where X is the vector superfield for the U(l) gauge supermultiplet and fei,2 3 are dimensionless 
couplings. There result corrections to the scalar potential: 

AVf = -^[(M+h)\Sm A + k 2 \S\ 2 \S\ 2 (\S\ 2 + \S\ 2 ) 2 }+... (3.26) 

AV ° = ^2(l 5 l 2 -^l 2 )^l 5 l 4 - fc 3^| 4 ) + --- ■ (3-27) 
The couplings k\, k 2 and run with RG scale according to 

(3(k 1 )=f3(k 3 ) = ^(3(k 2 ) = -^, (3.28) 



which, together with eqs. ( |3.26| ) and ( 3.27 ), explains the presence of the last terms in eq. ( 3.24 ). 



Note that far from the flat direction, the non-renormalizable D-term contributions in eq. (3.27) 
become comparable to all of the other terms in the scalar potential for intermediate-scale field 
values. 

Fortunately, one can see that at the minimum near the flat direction with \(S)\ = \{S)\ = 
s ~ yJmy/M, all of the "extra" terms in eq. ( |3.24| ) which were not reflected in eq. fl3.13|) are 
suppressed by at least one additional factor of order mw/M. In checking this, it is useful to 
note that the expectation value for the deviation from D-flatness is given by 

\{S)\ 2 - \{S)\ 2 = 7^(™ 2 ~ m 2 ) + ^(h ~ h)s\ (3.29) 

which is suppressed with respect to s 2 by at least mw /M. 

We have shown that the existence of an intermediate scale vacuum expectation value is 
determined by a one-loop effective potential which suffers only small corrections with respect 
to the tree-level value as long as renormalized couplings near the scale of the putative VEV are 
used. In particular, the conditions necessary for the existence of the intermediate scale VEV 
do not depend (to the lowest non-trivial order in mw/M) on non-quadratic Kahler potential 
parameters over which we have little control. On the other hand, eq. ( |3.29|) shows that the 
deviation from D-flatness, while appropriately suppressed, does depend on higher-order Kahler 
potential parameters which violate the approximate S <-> S symmetry; both terms in eq. (3.29) 



are of order rn^y. These D-term contributions are important since all scalars ifi in the low- 
energy theory obtain contributions to their squared masses equal to 

Aml % = qi g 2 (\(S)\ 2 -\(S)\ 2 ), (3.30) 
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where </j is the U(l) charge of </>j. This effect has been studied in numerous papers including 
[17, 18, 16, ^O], pl| , but the presence of the Kahler-potential-dependent terms in eq. dpi) does 



not seem to have been emphasized before. The presence of these terms makes it very difficult 
to predict the magnitude or even the sign of the D-terms in the low energy theory, since the 
Kahler potential terms are generally not constrained by symmetries. On the other hand, it also 
means that the existence of L>-terms in the low-energy theory should be quite generic, and a 
useful phenomenological prediction of the general scenario espoused in this paper. 

4 Minimal model of automatic imparity from gauged B — L 

One example of a model in which an hcfA term is used to generate an intermediate scale was 
proposed in |fl^| . The motivation for this model is as follows. The MSSM is usually defined 
to respect -R-parity (or equivalently, matter parity) which prevents rapid proton decay and 



provides for a stable neutralino LSP cold dark matter candidate [22|. However, in the MSSM 
the imposition of i?-parity is somewhat ad hoc, since failure to impose it would not result in 
any internal inconsistency. The conservation of matter parity (or equivalently, i?-parity) in the 
MSSM can be explained in terms of a deeper principle by starting with a continuous gauged 
B — L symmetry, where B is baryon number and L is the total lepton number. Since matter 
parity is defined to be 

P M = {-lf {B - L \ (4.1) 



a gauged U(1)b-l will forbid all -R-parity violating operators 24]. There is no massless 
gauge boson that couples to B — L, so this gauge symmetry must be spontaneously broken. 
However, provided that all VEVs or other order parameters in the theory carry 3(B — L) charges 
which are even integers, the Z2 matter parity subgroup must survive as an unbroken and exact 



remnant of the original gauged B — L symmetry [24, 25, 26, 27]. 

In a minimal extension of the MSSM that realizes this idea, one embeds U(1)y into U (1)r x 
U(1)b-Li introducing a single new gauge multiplet. To provide for the breaking 

U(1) R x U(1) B -l - U(1) Y x z 2 matterpaxity , (4.2) 

one also introduces two chiral superfields S, S with B — L charges —2 and +2 respectively. These 
fields play precisely the same role as their namesakes in section 3. In addition to breaking B — L 
to -R-parity, an intermediate-scale VEV for S provides a see-saw mechanism for neutrino masses 
by means of superpotential couplings 

W = \y l i SViUj + yl^HuLj . (4.3) 
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Here 17, (i = 1,2,3) are three gauge-singlet neutrinos; the ordinary Standard Model neutrinos 
V{ are contained in the S , f7(2)^-doublet fields Lj. (The Fj must exist anyway in order to provide 
for U(1)b-l anomaly cancellation.) The resulting 6x6 see-saw neutrino mass matrix in the 
(vi,Vi) basis is 



2 \VuYu (S)ys, 

where v u = (H%) = 174 GeV sin/3. In order to avoid a disastrously large U(1)b—L -D-term, the 
scalar component of S must get a VEV which is very nearly equal to (S). Many prominent 
features of this model, including the mass spectrum of the S, S supermultiplets, have already 
been studied in some detail in ref. [14], Here, I will briefly review the results found there, and 
extend the analysis by examining the solutions to the RG equations for the /i-couplings to 
explore which parts of parameter space can actually lead to symmetry breaking with (S) ~ (S) 
at an intermediate scale. 

The superpotential, supersymmetric scalar potential, and supersymmetry-breaking terms 
are all just given by the same expressions eqs. ([O]), (p^ ) and (|3.3|) as in section || (with 
g — > gx given below). In fact, the important features of this model are all identical to the ones 
already discussed there, except that the B — L gauge charges are normalized differently, and 
there is a non-trivial kinetic mixing due to the presence of two U(l) gauge symmetries. In the 
U(1)y, U(1)b-l basis, S and S have charges (0, ±2). Another useful basis is the choice U(1)r, 
U(1)b-l, under which S and 5 have charges (1,-2) and (—1,2). However, neither of these 
bases avoids kinetic mixing in this model. In general, the couplings of the ?7(1) gauge bosons to 
a matter field fa with U(1)r and U(1)b—l charges Ri and (B — L)i are given by the covariant 
derivative 

Dpfa = (d^ + ig?AK + igf- L A B - L )fa, (4.5) 



9f 
9? 



g R Ri + gB~L,R\f3/8(B-L)i, (4.6) 
g B -L^/8(B - L)i + g R , B -LRi- (4.7) 



where gB-L, 9r, 9b-l,r and gR^B—L are distinct gauge couplings. In terms of these couplings, 
the weak hypercharge coupling gy and the coupling gx which appears in the D-term scalar 
potential V D = \g 2 x {\S\ 2 - \S\ 2 ) 2 are given by 

9x = (9r-\/V29b-l,r) 2 + (^3/29b-l-9r,b-l) 2 , (4.8) 
g Y = \f5/2(gRgB-L - 9b-l,r9r,b-l)/9x- (4.9) 
If one starts at the unification scale with S'O(10)-like boundary conditions 

93 = 92 = 9R = 9B-L = 9V, 9B-L,R = 9RB-L = 0. (4.10) 
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near Qq ~ 3 x 10 16 GeV, then the successful low-energy prediction for sin 2 9\y is maintained 
after the symmetry breaking In that case one can most easily run RG equations by going to 
the one- loop multiplicatively renormalized basis, as shown in ref. [14]. 

For the boundary conditions eq. ( 4.1C| ), the analytic one-loop results for quantities of interest 
are as follows. The running gauge couplings at any scale t are given by 



9R = 9B-L = 9u(k+ + «-)/2; 
gn,B-L = 9b-l,r = 9u(k- - K+)/2, 



(4.11) 
(4.12) 



where 



K± 



l + |^(9±V6)(*o-t) 



-1/2 



(4.13) 



are monotonically decreasing as t decreases into the infrared. [The first equality in each of 
eqs. ( PH) and dO^ ) is due to an amusing but otherwise inessential coincidence in the matter 
content of this model.] Now let us assume for simplicity that near the gauge coupling unification 
scale to, the supersymmetry breaking masses obey supergravity inspired boundary conditions, 
with a common scalar ttjq, a common gaugino mass m.i/2> an d a specified Hq/xq. The running 
gaugino mass matrix in the R, (B — L) basis is given by 



M, 



1 

2 m V2 



(4.14) 



If the Yukawa couplings in eq. (4.3) are small enough to be neglected in the RG running, then 
the squared masses for the scalar components of S, S run according to 



2 2,2 



'33 + W6\ (1 _ 4 



) + 



33 



(1 - K > 

150 J y 



The coupling x runs according to 



-(66+26-v/6)/75 -(66-26\/6)/75 
Xq K, K_ , 



(4.15) 



(4.16) 



and, more importantly, the ratio h/x runs according to 



h/x = h /x - (t Q - t)|^mi /2 [(5 + 2y/6)n% + (5 - 2y/E) 



(4.17) 

In order to decide whether a non-trivial intermediate scale minimum (S) ~ (S) will occur, 
one should compare the running values of \h/x\ 2 and m| + m~ at a scale Q ~ (5) ~ v /m vy-^P ) 
as described in section ||. The running of \h/x\ 2 as a function of scale is shown in Figure [2|, 
with a factor of w- 2 ^ 2 removed to make it dimensionless. For illustrative purposes, I have chosen 
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3 4 5 
Log 10 [Q /Q] 



Figure 2: RG running of h/xm\j2 for the minimal model of automatic -R-parity conservation 
described in the text. Initial values at the high scale, corresponding to the left-hand side of the 
graph, are ho/xomi/2 = 0, —1, —2, —3 (solid lines) and 1, 2, 3 (dashed lines). 

initial values /io/ 3; o m i/2 = 0] if] ±3 at the unification scale, with g^/An = 0.04. The fi gure 
illustrates that for negative values of ho/xQirii^, there is constructive interference that allows 
\h/x\ 2 to grow in the infrared (to the right of the figure). The rate of running gradually slows 
down as the U{\) gauge couplings diminish at lower scales. Conversely, positive initial values 
for /io/xo ?Tl i/2 l ea d to a destructive interference in the magnitude of \h/x\ 2 at lower scales. 

Of course, the scalar masses m 2 s and w| are also running because of the £7(1) gauge cou- 
plings. This effect is quite significant because of the large (±2) values of B — L charges of the the 
S, S fields, and it opposes the symmetry breaking. Therefore, there is a competition between 
the magnitude of running for m| + m|- and \h/x\ 2 to decide whether an intermediate scale VEV 
can occur. As in section 3, this requires \h/x\ 2 > 6(m| + mJ|) at an intermediate scale. For a 
numerical example, one can use gfj/^ir = 0.04 and Qo/Q = 10 6 ov to — t = 13.8 in eqs. Q4.15 ) 



and ( 4.17] ) to get m 2 s = m 2 ^ « tuq + 0.33m^ 2 and h/x 



ho/xo — Q.&Smi^- In Figure ^ I show 
the regions of parameter space in the /io/xow-i/2 vs. ?7ig/m^ 2 plane for which a non-trivial 
local minimum occurs along the S, S flat direction in this model. There is a local minimum for 
all values of m^/m 2 ^ below the lines indicated. The solid line corresponds to a ratio of the 
high scale to the intermediate scale of 10 6 , while the dashed line corresponds to the case of no 
renormalization of parameters (if the intermediate scale were pushed up close to the high scale, 
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-5 -4 -3 -2 -1 1 2 3 4 5 

h /x m 1/2 

Figure 3: The parameter space in the /io/xo m i/2 vs. mQ/fflL plane for which a non-trivial 
local minimum occurs along the S, S flat direction in the model described in the text, neglecting 
Yukawa couplings. A local minimum exists for all values of rog/m^ below the lines shown. The 
solid line corresponds to a hierarchy of the high scale to the intermediate scale of Qo/Q = 10 6 , 
while the dashed line corresponds to the case of no renormalization (Qo/Q = 1)- 

which can occur in the limit that x is very small). As Fig. |3| shows, the parameter space in 
which intermediate-scale symmetry breaking occurs is actually diminished by RG running for 
ho/xomi/2 ^ — 2, and is increased for /io/xo m i/2 ~ — 2. If one assumes a boundary condition 
m i/2 = ^ m o as can happen in certain dilaton-dominated scenarios of supersymmetry breaking, 
then a non-trivial minimum at an intermediate scale requires /io/£o m i/2 ~ — 1-9 or > 3.7 in 
this picture. 

The discussion above is somewhat conservative, in that we have not included possible neg- 
ative RG contributions to m| and w|. One might imagine that one could increase the region 
of parameter space in which symmetry breaking occurs by including the effects of the Yukawa 
couplings y v and ys in eq. ( [4.3[ ). The latter coupling will certainly act to decrease m 2 s in the 
infrared, which clearly favors a VEV for S. However, the impact of these corrections is some- 
what limited. For example, assuming one neutrino Yukawa coupling ys and its corresponding 
scalar cubic term as dominate, the relevant beta functions are: 

WTr 2 (3(m s ) = \ys\ 2 (m s + 2rr4) + \as\ 2 -8j29 2 a MlC a (S) (4.18) 

a 
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16vr 2 /3(m|) = 2\y s \ 2 (m 2 s + 2m|) + 2|a s | 2 + . . . - 8]T g 2 a M 2 C a (v) (4.19) 

a 

where the ellipses in eq. Q4.19Q refer to additional positive-definite contributions from Yukawa 
couplings. A comparison of the terms in these two equations reveals that without additional 
interactions, it is quite difficult to arrange for m| to be driven negative without first driving 
negative, since C a (S) = 4C a (F). Once this occurs, m 2 s will receive large positive contributions 
from the 2\ys\ 2 rn 2 7 term in its beta function. One possible alternative is to allow S and S to have 
superpotential Yukawa couplings to some heavy particles which also have gauge interactions 
with respect to some strongly-coupled gauge group. The soft squared masses of these new 
heavy scalars will be large and positive, and through Yukawa couplings can allow m 2 s and/or 
m 2 ^ to be driven negative or at least smaller in the infrared, thus increasing the parameter space 
in which intermediate-scale symmetry breaking occurs. In any case, the contributions of the 
dimensionless super symmetry-breaking terms are likely to be non-negligible, and always favor 
symmetry breaking. 

5 Model with automatic imparity conservation tied to a solu- 
tion to the fi problem 

In this section I consider an extension of the model in the previous section which also incor- 
porates a solution to the /z-problem together with an invisible axion of the type proposed in 
ref. |2S} |. In addition to the U(1)b-l gauge supermultiplet and S,S fields, introduce an addi- 
tional pair of neutral chiral supermultiplets N and N. These fields are charged under a global 



anomalous Peccei-Quinn (PQ) symmetry [29| as shown in Table 2. When the fields (N) and 
(N) get VEVs, they will spontaneously break this PQ symmetry, which is also explicitly bro- 
ken by the QCD anomaly. As before, I will assume that superpotential mass terms are absent. 



Table 2: Charges of Higgs fields in the model described in section [|. The B — L and Y 
symmetries are local, and the PQ symmetry is a global symmetry with a QCD anomaly. 





S 


S 


N 


N 


H u 


H d 


B-L 


2 


-2 














Y 














1/2 


-1/2 


PQ 


1 


-1 


-1 


1 


1 


1 
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The leading terms in the superpotential consistent with the symmetries in Table 2 can then be 
written as: 

W = ^jS 2 S 2 - ^pSSNN + ^jN 2 N 2 + ^N 2 H u H d . (5.1) 
2M M 2M M y J 



The last term in eq. fl5.1|) becomes the term of the MSSM when iV obtains its VEV. Then 



(N) should be of order the intermediate scale yJmy/M for two distinct reasons; first in order to 
allow the effective contribution to \i = x fl (N) 2 /M to be of order mw, and second to allow the 
PQ scale to fall within the window permitted by direct searches and astrophysical constraints 

Because we are looking for a minimum of the potential with (S), (S), (N), (N) 3> (H u ), {Hd), 



the last superpotential term in eq. (5.1) is only a small perturbation on the dynamics that fixes 
the VEVs of S, S, N, and N. The relevant supersymmetric part of the scalar potential is 
therefore 

^SUSY = ^(\S\ 2 + \^ 2 )\ Xl SS-x 2 NN\ 2 + ^m 2 + \N\ 2 )\x 3 NN-x 2 S^ 2 

+ 9 f(\S\ 2 -\S\ 2 ) 2 , (5-2) 
and the supersymmetry-breaking terms are 

^breaking = ( — S 2 ^ 2 - — SSNN + ^N 2 N 2 \ + C.C. 
breaking \2M M 2M J 

+m 2 s \S\ 2 + m||5| 2 + m^|iV| 2 + m^N^ . (5.3) 

The parameter space of this model is multidimensional and complicated, so I will be content 
to explore it quantitatively in some special cases. In order to simplify the analysis, I will assume 
that x 2 jx\ and x$/x\ are both real. Then each of x\, x 2 , and X3 can be made real and positive 
by a field rephasing, without loss of generality. I will also assume that 

h\jx\ = h 2 /x 2 = h 3 /x 3 = (h/x) (5.4) 

at the unification scale Qo, and that all scalars have a common squared mass equal to m 2 , at that 
same scale. RG equations are used to run the parameters hi/ximi/2, h 2 /x 2 rrii/ 2 , m|/m 2 ^ 2 , and 
m^/rriy 2 to the intermediate scale Q to evaluate the scalar potential. The ratios h 3 /x 3 mi/ 2 , 
m 2 N /m 2 y 2 , and m jf/ m i/ 2 involve gauge singlets and do not run. Then minima of the potential 
can be identified by a straightforward procedure. This involves numerically solving quartic and 
quadratic equations to find all possible candidate extrema, and then requiring that all of the 
eigenvalues of the Hessian matrix at the minimum are positive. The candidate minima occur 
near the _D-flat direction \(S}\ = \(S)\ = s. Furthermore, my assumption m 2 N = m 2 ^- implies 
that \{N)\ = \ {N)\ = n, because of an approximate symmetry N N. 
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I have checked that for generic choices of xi/x% and X2/X3, there are stable local minima 
for sufficiently large values of {h/x)o/m\/2 compared to m§/fflw 2 , just as for the model in the 
previous section. However, in the present case both s and n are non-zero. It is instructive to 
first consider the special case that Qo/Q = 1 so that there is no renormalization of parameters 
from the high scale. Then, at least for many values of X1/X3 and X2/X3, a stable local minimum 
exists provided that 

{h/x)l > 12ml, ( 5 - 5 ) 

with n 2 /s 2 = (x\ + X2)/(x3 + X2)- For the more realistic case including renormalization effects, 
I have chosen to present results for the special case xi/x% = 4 and X2/X3 = 3, and for case 
xi/%3 = 1/4 and X2/X3 = 1/3, in Fig. (These ratios are taken to be specified at the 
intermediate scale.) For both models, the area below the dashed curve is the region of parameter 
space in which a local minimum exists when Qo/Q = 1, as in eq. Q5,5| ). The area below the dash- 
dot line is the corresponding region for Qo/Q = 10 6 when (xi/x^, X2/X3) = (4,3). Likewise, the 
area below the solid line is the region allowing intermediate scale VEVs for Qo/Q = 10 6 and 
(X1/X3, X2/X3) = (1/4,1/3). I have checked that similar results obtain for many other choices 
of X1/X3 and X2/X3. 

As Fig. H shows, the critical line for (/i/x)o/ m i/2 < is surprisingly insensitive to renor- 
malization group running. This is actually something of an accident, reflecting the fact that 
the scalar squared masses and (hi/xi) 2 run at comparable rates. For the opposite sign of 
{h/x)o/mi/2, there is somewhat more sensitivity. Also, I find that for values of m 2 , far below 
the critical lines shown, there are sometimes several distinct non-trivial local minima separated 
by "ridges" in field space. Note that the critical lines in this model are between the two lines 
shown in Figure ||. This can be understood as follows: along the flat direction (n, s) chosen by 
the model in this section, the minimum is determined by effective h and x couplings which are, 
roughly speaking, "averages" of some couplings that are renormalized by B — L interactions just 
as in section 4, and other couplings that are not renormalized because they are gauge singlets. 

However, it is important to realize that the situation depicted in Fig. |] is not inevitable, since 

the relationship between hi, hi, and /13 can be relaxed. To see this, note that if the special 

(one-loop RG-invariant) relationship X1X3 = x\ were to hold, then there would be exactly 

flat directions^ of the supersymmetric part of the potential, eq. ( |5.2[ ), with |5| 2 = |5| 2 and 

NN/SS = x\/x2 = X2/X3. Now, with the boundary conditions on the hi/x{ discussed above, it 

happens that at one loop order the supersymmetry-breaking quartic part of the scalar potential 

is also flat at one-loop order. However, if one imagines relaxing this condition and allowing 

^Of course, these flat directions will be lifted by higher-order terms in the superpotential, but those terms are 
suppressed by an additional factor of 1/M 2 . 
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Figure 4: The parameter space in the /io/^o m -i/2 y s- mo/ m i/2 pl ane f° r which a local minimum 
involving non-zero (S) = (S) and (iV) = (N) occurs, for the model and RG boundary conditions 
described in section |5|. At least one non-trivial local minimum exists for all values of m,Q/m^ 2 
below the lines shown. For a hierarchy of the high scale to the intermediate scale of Qo/Q = 10 6 , 
the dash-dot line corresponds to the choice X1/X3 = 4 and X2/X3 = 3 and the solid line 
corresponds to x\jx^ = 1/4 and X2/X3 = 1/3. In the case of no renormalization (Qo/Q = 1), 
both models have non-trivial local minima for all m^/m^ below the dashed line. 

generic values of hi/xi, then the negative supersymmetry-breaking /i-terms will dominate the 
scalar potential for field values slightly above the intermediate scale for some neighborhood of 
parameter space near X1X3 = x\. This makes it clear that for even for rather small, but generic, 
values of (h%, /12, /13), there must always be at least some finite region of (x±, X2, X3) parameter 
space near X1X3 = x\ in which the symmetry breaking does take place. I have checked that 
this region of parameter space can be quite large. 

The model outlined above has the nice feature that it relates the PQ intermediate scale 
governed by (N) to the neutrino see-saw intermediate mass scale governed by (S). A similar 
model was proposed in ref. |3l|], but in that paper the mechanism for intermediate-scale sym- 
metry breaking was supposed to be negative RG corrections to scalar squared masses. In the 
model described here, B — L is gauged in order to provide for automatic imparity conservation, 
and this makes it extremely difficult for negative scalar squared masses m 2 s and to arise. 
In any case, I have argued that the hcj) 4 coupling probably provides a significant part of the 



23 



effect. If the physical neutrino masses are very small, it may be required that there is a mild 
hierarchy (N) < 10 12 GeV <C (S) in order to accommodate the allowed axion window. This 
can be accomplished by choosing the couplings x\,X2, x% appropriately. In the model I have de- 
scribed, the intermediate scale symmetry breaking occurs even though no scalar squared mass 
runs negative, and without tuning the relative magnitudes of any superpotential terms. As 
is the case for the model in the previous section, matter parity is automatically a conserved 
symmetry because of the way B — L is broken. In addition, the strong CP problem is solved 
and there is an invisible axion (which is mainly lm[N — N] / y/2) along the lines of ref. 28]. The 
fermionic axino and scalar saxino partners of the axion both obtain electroweak-scale masses, 
but they are only very weakly coupled to ordinary matter. As before, one expects that the 
pattern of MSSM squark and slepton masses will be augmented by a B — L D-term contribu- 
tion, which could be distinguished at the CERN Large Hadron Collider or a future lepton linear 
collider. As emphasized in section ||, the magnitude and even the sign of the D-term depend 
on non-renormalizable Kahler potential terms, but its existence is non-negotiable. 



6 Conclusions 

In this paper, I have studied the effects of super symmetry-breaking <ft 4 couplings in producing 
intermediate-scale VEVs. Any truly complete model of supersymmetry breaking should be able 
to predict the magnitude of these terms, at least relative to the corresponding superpotential 
couplings. As shown in sections 2 and 3, these terms are most important in or near D-flat 
directions that are also F-Qat at the renormalizable level. For these purposes, there is a well- 
defined procedure for including RG effects, as given by the beta function eq. ( |2.8| ), which 
incorporates the large logarithms in an effective potential approach. 

An intermediate scale roughly of order V m VK -^Planck is suggested both by the see-saw sce- 
nario for neutrino masses and the allowed axion window (3^]. In sections 4 and 5, I discussed 
simple models in which the existence of these scales is dependent on the presence of dimen- 
sionless supersymmetry-breaking terms. The numerical studies in this paper show that in the 
simplest cases, one requires the values of the dimensionless ratio h/xm^/v at the Planck scale 
to be fairly large, and perhaps to favor a particular sign, in order to achieve intermediate-scale 
symmetry breaking. However, this may well be too conservative since Figs. || and |] neglect 
the possible effects of negative contributions to scalar squared masses from RG effects, even 
if the scalar squared masses never run negative. Furthermore, as I argued briefly in section 
5, the parameter space in which intermediate-scale symmetry breaking can take place can be 
increased substantially if one is close to an accidental nearly-flat direction and the h(f> terms 
are misaligned with the corresponding superpotential terms. In any case, it is likely that the 
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dimensionless super symmetry breaking couplings play a crucial role in producing VEVs along 
supersymmetric flat directions. 
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Appendix: flat directions in the MSSM 

In this paper, I have mainly been concerned with the ability of dimensionless supersymmetry- 
breaking couplings to produce desired intermediate-scale VEVs like the PQ scale and the neu- 
trino seesaw scale. It is also interesting to consider the opposite case, when we do not want 
an intermediate scale VEV to occur; namely, /i-terms corresponding to flat directions in the 
MSSM. Supersymmetric flat directions are parameterized by analytic gauge-invariant polyno- 
mials in chiral superfields. For example, there are flat directions QQQL and Wide which are not 
lifted by any renormalizable supersymmetric terms in the scalar potential. Instead, they are 
lifted by super symmetry-breaking (mass) 2 terms, by non-renormalizable superpotential terms, 
and by the corresponding h-type terms. The latter terms always lower the scalar potential for 
some choice of phases of the fields, and therefore favor the existence of a non-trivial minimum 
at an intermediate scale. Naively, these minima are potentially dangerous, since they of course 
break color and electric charge. However, the mere existence of such vacua does not imply that 
the universe has to be in one of them in the present day, even if they are global minima. On 
the other hand, the dynamics associated with the flat directions may play an important role in 
the early universe, particularly during inflation when the relevant F-term may be the auxiliary 
component of the inflaton, and in the study of baryogenesis [32, 11]. Here, I will illustrate the 



importance of the RG running by considering the /i-terms for the uude flat directions. The 
analysis is very similar e.g. for the QQQL flat directions. 
The relevant superpotential is given by 

^ijkl _ 

W = -^UiUjd k ei (A.l) 

where gauge indices are suppressed in the only possible way and i = 1,2,3 is a family index. 
The super symmetry-breaking potential is given by 
foijkl '> 



V 



UiUjd k ei + c.c. + ( m «J u i| 2 + m2 d r \di\ 2 + m^M 2 ) • ( A - 2 ) 



2M 

i=i 



Now, suppose we are interested in a particular flat direction with SU(3)c color components 

'ON _ /0 

u x = | | ; u 2 =\<j>\; di = I | ; ei = 0. (A.3) 
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Here the indices 1 and 2 on u\, U2, d\ and e\ represent arbitrarily-labelled family components. 
The potential along the flat direction is given by 

V = E-Vl 2 + (A^4 + c . c .) + (A .4) 

where 



^m 2 = m 2 Ui + m 2 U2 + m 2 di + m 2 ei ; (A.5) 



x 2 



h 1 ' 11 ; (A.6) 

I (V 211 | 2 + |x im | 2 + \x 12a \ 2 + \x 12U \ 2 ) . (A.7) 
i=l 

In the MSSM, there are no Yukawa couplings which simultaneously link any pair of fields u, 
d and e. Furthermore, to a very good approximation, the Yukawa contribution to wavefunction 
renormalizations are only non-vanishing for the t, b, r family. Therefore, the conditions leading 



to eq. (|2.10|) are satisfied, and one can write a simple RG equation for each separate ratio 

>^) = ^ + ^M, (A.8) 

(I emphasize that this is true even when the coupling involves t, b, r fields.) The resulting run- 
ning of this family-independent ratio is graphed in Fig. [|, starting from unified gauge couplings 
at Qo = 3 x 10 16 GeV. As the graph shows, the ratio {h/xnii^) 2 can become quite large for 
negative ho/xomi/2, while for positive ho/xomi/2 it tends to run to very small values. 

A non-trivial local minimum in the scalar potential can occur at an intermediate scale if 

\h\ 2 -3x 2 J]m 2 > 0. (A.9) 

Note that while the RG running shown involves only the coupling x corresponding to the 
gauge-invariant polynomial of the flat direction, the minimization condition that needs to be 
satisfied involves the larger averaged squared coupling x 2 . If only one coupling dominates, then 
x 2 ~ x 2 . The relationship between x 2 and x 2 is quite model-dependent; however, there are 
strong constraints on these couplings because they are dimension 5 proton decay operators. 
Therefore, it is probable that the largest of the couplings x belong to the t, b, t family, which 
are constrained the least by proton decay searches. In addition, the top and bottom squarks 
get negative RG corrections to their masses due to the top and bottom Yukawa couplings. 
Therefore, it is quite likely that if a non-trivial local minimum exists, it will involve the third- 
family quarks and sleptons. I find that such a minimum can easily exist for h/xm.1/2 ^ — 2 or 
so, because of the fast growth of h/x in the infrared. However, the precise regions of parameter 
space in which this can occur depend on numerous unknown variables (including the scalar 
cubic couplings) so I will not attempt to outline them here. 
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Figure 5: RG running of h/xm\/2 for a wade fiat direction in the MSSM as described in the 
text. Initial values at the high scale (corresponding to the left-hand side of the graph) are 
h/xmi/2 = 0,-1,-2,-3 (solid lines) and h/xmi/2 = 1,2,3 (dashed lines). 

References 

[1] H.P. Nilles, Phys. Rept. 110, 1 (1984); H.E. Haber and G.L. Kane, Phys. Rept. 117, 75 
(1985). 



[2] 
[3] 
[4] 
[5] 

[6] 
[7] 
[8] 

[9] 
[10] 



S.P. Martin, "A Super symmetry primer, " |hep-ph /9709356; . 
K. Harada and N. Sakai, Prog. Theor. Phys. 67, 1877 (1982). 
L. Girardello and M.T. Grisaru, Nucl. Phys. B194, 65 (1982). 

K. Inoue, A. Kakuto, H. Komatsu and S. Takeshita, Prog. Theor. Phys. 68, 927 (1982), 
erratum ibid 70 330 (1983). 

D.R.T. Jones, L. Mezincescu and Y.P. Yao, Phys. Lett. 148B, 317 (1984). 
L.J. Hall and L. Randall, Phys. Rev. Lett. 65, 2939 (1990). 

F. Borzumati, G.R. Farrar, N. Polonsky and S. Thomas, "Soft Yukawa couplings in super- 
symmetric theories," hep-ph / 9902443| . 



I. Jack and D.R.T. Jones, "Nonstandard soft supersymmetry breaking," bep-ph /9903365| . 



F. Buccella, J.P. Derendinger, S. Ferrara and C.A. Savoy, Phys. Lett. 115B, 375 (1982); 
R. Gatto and G. Sartori, Phys. Lett. 157B, 389 (1985); M.A. Luty and W.I. Taylor, Phys. 
Rev. D53, 3399 (1996); 



27 



[11] M. Dine, L. Randall and S. Thomas, Nucl. Phys. B458, 291 (1996). 

[12] T. Gherghetta, C. Kolda and S.P. Martin, Nucl. Phys. B468, 37 (1996). 

[13] G. Cleaver et al, Nucl. Phys. B545, 47 (1999). 

[14] S.P. Martin, Phys. Rev. D54, 2340 (1996). 

[15] G. Lazarides and Q. Shan, Phys. Rev. D58, 071702 (1998). 

[16] I. Jack et al, Phys. Rev. D50, 5481 (1994). 

[17] M. Drees, Phys. Lett. 181B, 279 (1986). 

[18] J.S. Hagelin and S. Kelley, Nucl. Phys. B342, 95 (1990); A.E. Faraggi, J.S. Hagelin, 
S. Kelley and D.V. Nanopoulos, Phys. Rev. D45, 3272 (1992). 

[19] Y. Kawamura and M. Tanaka, Prog. Theor. Phys. 91, 949 (1994); Y. Kawamura, H. Mu- 
rayama and M. Yamaguchi, Phys. Lett. B324, 52 (1994); Phys. Rev. D51, 1337 (1995). 

[20] H.C. Cheng and L.J. Hall, Phys. Rev. D51, 5289 (1995). 

[21] C. Kolda and S.P. Martin, Phys. Rev. D53, 3871 (1996). 

[22] G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996); J. Ellis, 
T. Falk, K.A. Olive and M. Schmitt, Phys. Lett. B413, 355 (1997); J.D. Wells, "Mass 



density of neutralino dark matter," hep-ph /9708285 , in "Perspectives on Supersymmetry" , 
ed. G.L.Kane, World Scientific. 

[23] R.N. Mohapatra, Phys. Rev. D34, 3457 (1986); A. Font, L.E. Ibanez and F. Quevedo, 
Phys. Lett. B228, 79 (1989). 

[24] S.P. Martin, Phys. Rev. D46, 2769 (1992). 

[25] R. Kuchimanchi and R.N. Mohapatra, Phys. Rev. Lett. 75, 3989 (1995); Z. Chacko and 
R.N. Mohapatra, Phys. Rev. D58, 015001 (1998). 

[26] K. Huitu, P.N. Pandita and K. Puolamaki, Phys. Lett. B415, 156 (1997). 

[27] C.S. Aulakh, A. Melfo and G. Senjanovic, Phys. Rev. D57, 4174 (1998); C.S. Aulakh, 



A. Melfo, A. Rasin and G. Senjanovic, |hep-ph/9902409 



[28] J.E. Kim and H.P. Nilles, Phys. Lett. 138B, 150 (1984). Mod. Phys. Lett. A9, 3575 (1994). 
[29] R.D. Peccei and H.R. Quinn, Phys. Rev. Lett. 38, 1440 (1977). 

[30] J. Preskill, M.B. Wise and F. Wilczek, Phys. Lett. 120B, 127 (1983); L.F. Abbott and 
P. Sikivie, Phys. Lett. 120B, 133 (1983); M. Dine and W. Fischler, Phys. Lett. 120B, 137 
(1983); M.S. Turner, Phys. Rept. 197, 67 (1990); G.G. Raffelt, Phys. Rept. 198, 1 (1990). 

[31] H. Murayama, H. Suzuki and T. Yanagida, Phys. Lett. B291, 418 (1992). 

[32] I. Affleck and M. Dine, Nucl. Phys. B249, 361 (1985). 



28 



